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i 

SUMMARY 

I 

k  An  analytical  and  experimental  study  was  conducted  to 

*  determine  criticality  of  interlaminar  disbonds  by  NOG  methods. 

Criticality  of  such  flaws  in  a  shear  environment  (action  of 
shear  near  support)  was  defined  in  terms  of  crack  propaga¬ 
tion  and  was  analyzed  by  principles  and  methods  of  fracture 
mechanics.  Growth  of  disbonds  under  cyclic  loading  was  also 
studied.  Failure  under  compressive  loading  in  presence  of 
a  disbond  was  defined  in  terms  of  buckling  and  an  elastic 
stability  analysis  was  utilized  for  assessing  criticality. 

Analytical  predictions  were  compared  with  experimental  results 
in  both  cases.  Ultrasonic  "C"  scams  were  used  for  detection 
and  tracking  of  the  flaws.  Preliminary  wave  propagation 
|>  studies  were  conducted  for  estimating  changes  in  storage  ] 

and  loss  moduli  induced  by  moisture  conditioning  to  examine  ; 

the  possibility  of  using  such  measurements  as  NDT  techniques 
I  for  damage  assessment. 
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INTRODUCTION 

The  NDE  analyst  faces  problems  and  challenges  in  evaluating 
flaws  in  composite  structures  because  of  the  complexity  of  fiber 
composites.  At  the  present  time  only  limited  use  is  made  of  NDT 
techniques  in  quantitative  evaluation  of  life-lin,  ‘■-ing  composite 
material  flaws.  This  is  due  to  limitations  both  the  state  of 
the  art  in  assessing  flaw  criticality  in  composite  structures  and 
of  contemporary  NDE  techniques.  However,  as  the  understanding  of 
flaw  criticality  in  composites  and  new  NDE  techniques  are  evolved, 
more  reliable  serviceability  criteria  may  be  adopted.  This  will 
enable  definition  of  regular  inspection  intervals  for  a  structural 
element  during  its  service  to  determine  the  criticality  of  exist* 
ing  flaws. 

Development  of  the  required  nondestructive  evaluation  method* 
ology  includes  consideration  of  the  following  successive  phases  of 
interests 

1.  Definition  of  measurable  quantities.  This  defines  what 
can  be  measured  at  the  lamina  and/or  laminate  level 
through  available  techniques. 

2.  Definition  and  demonstration  of  the  correspondence  be* 
tween  the  NDE  measured  quantities  and  the  associated  de¬ 
fects.  In  general,  the  defects  may  be  birth  defects  or 
service  defects,  and  the  effect  of  each  type  on  the  meas¬ 
urable  quantities  should  be  determined,  both  analytically 
and  experimentally. 

3.  Evaluation  of  the  criticality  of  the  defects.  This  re¬ 
quires:  evaluation  of  the  appropriate  residual  property 
(stiffness,  strength,  etc.)  in  the  presence  of  known 
defects;  an  experiment-analysis  correlation  study  for  the 
same;  and  formulation  of  a  quantitative  relationship  be¬ 
tween  the  NDE  measurements  and  these  residual  properties. 

In  the  present  program,  development  of  this  evaluation  method¬ 
ology  was  undertaken  for  defects  in  the  form  of  interlaminar  dis¬ 
bonds.  This  defect  was  selected  after  consideration  of  a  range  of 
possible  defects.  Various  kinds  of  defects  can  cause  relevant 
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strength  and  stiffness  degradation  in  composite  laminates;  namely; 
interlaminar  disbonds;  cracks  or  other  through- the- thickness  de¬ 
fects;  defects  in  bolted  and  bonded  joints;  damage  resulting  from 
impact;  and  fatigue  damage.  Different  NDI  techniques  (active  and 
passive)  which  are  usually  used,  or  being  developed  for  studying 
effects  of  such  defects,  are  listed  below; 

1.  Ultrasonics  (modulus  degradation  measurement  and  damage, 
detection) . 

2.  Acoustic  emission  (sequential  recording  of  the  damage 
growth  process) . 

3.  X-ray  and  thermography  (visual  and  real-time  detection 
of  defect  growth) . 

4.  Structural  vibrations  (stiffness  degradation  measurement) . 

5.  Penetrants. 

6.  Holography. 

For  the  present  work,  the  decision  was  made  to  use  ultrasonic 
"C"  scans  for  detecting  and  following  the  growth  of  interlaminar 
disbond.  This  defect  is  a  realistic  one  which  is  well-defined 
and  can  have  an  important  effect  upon  structural  performance.  The 
NDE  technique  is  widely  used  and  provides  a  reliable  measure  of 
the  size  of  the  interlaminar  disbond. 

Analytical  studies  were  conducted  to  model  the  physically  real¬ 
istic  interlaminar  damage  modes  and  to  assess  the  resulting  property 
degradation  and  criticality  of  disbonds  in  both  shear  and  compres¬ 
sion  environments.  Mechanical  tests  were  conducted  to  relate  the 
observed  states  of  magnitude  and  geometry  of  damage  to  the  residual 
performance  capability. 

In  addition,  preliminary  wave  propagation  studies  were  per¬ 
formed  for  quantifying  changes  in  storage  and  loss  moduli  Induced 
by  fatigue  damage  and/or  moisture  conditioning  and  to  examine  the  use¬ 
fulness  of  such  measurements  as  NDT  techniques  for  assessing  quan¬ 
titatively  the  magnitude  of  such  damage. 

Correlations  of  data  from  NDI,  mechanical  tests,  and  analyses 
were  carried  out  to  demonstrate  the  feasibility  of  using  this  ap¬ 
proach  for  the  development  of  quantitative  NDE  of  flaw  severity  in 
composites. 
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SCOPE 

Initiation  .and  growth  of  flaws  in  a  composite  structural  mem¬ 
ber  depend  on  imposed  loads  and  various  material  and  geometric  pa¬ 
rameters.  Applied  loads  can  be  classified  in  various  categories, 
namely : 

1.  Inplane  loads  (tensile,  compressive ,  shear  or  multi- 
axial)  ; 

2.  Transverse  loads  (and  shear) ; 

3.  Stresses  induced  by  environmental  factors  like  tempera¬ 
ture  and  moisture. 

Some  of  the  important  material  and  geometric  parameters  are: 

1.  Anisotropic  material  properties  and  orientation  of  in¬ 
dividual  laminae; 

2.  Laminate  layup; 

3.  Geometry  of  the  laminated  components;  and 

4.  Location,  shape  and  size  of  the  flaw. 

Also,  different  types  of  flaws  may  result  in  a  large  number  of  po¬ 
tential  failure  modes.  Therefore,  the  factors  which  influence 
the  phenomena  of  flaw  growth  identify  a  problem  of  enormous  com¬ 
plexity.  The  objective  of  the  present  program  was  to  choose  some 
representative  flaws  which  are  of  crucial  importance  in  composite 
structural  members  and  to  develop  a  methodology  for  quantifying 
damage  tolerance  characteristics  of  such  flawed  members  subjected 
to  the  type  of  loads  which  are  likely  to  cause  damage  growth  and 
final  failure.  Damage  tolerance  can  be  quantified  in  terms  of 
residual  stiffness,  strength  or  lifetime. 

The  procedure  employed  was  to  combine  analytical  methods  for 
assessing  criticality  of  flaws  and  residual  properties  with  an  ex¬ 
perimental  program  which  supplemented  the  analytical  investigation 
The  other  alternative  was  an  extensive  experimental  investigation 
which  could  become  cost  prohibitive  even  for  a  single  type  of  flaw 
It  was  also  noted  that  it  would  be  impossible  to  use  a  single  ana¬ 
lytical  model  to  encompass  the  whole  spectrum  of  fatigue,  fracture 
or  other  kinds  of  failure  for  various  types  of  flaws.  A  logical 
first  step  was  to  begin  with  a  simple  but  commonly  encountered 


NADC-77278-30 


flaw  geometry  like  an  interlaminar  disbond  in  a  laminated  beam  and 
to  carry  out  analytical  and  experimental  correlation  studies,  so 
that  confidence  could  be  developed  for  similar  studies  on  other 
types  of  flaws.  Interlaminar  disbonds  are  quite  common  in  lamina 
ated  members  because  they  can  exist  as  "birth"  defects  and  can  be 
created  during  service  by  foreign  object  damage  or  various  other 
reasons.  Laminated  members  are  used  for  carrying  inplane  loads  as 
well  as  transverse  shear  stresses.  Presence  of  a  disbond  does  not 
affect  the  primary  function  of  a  member  under  inplane  loads  except 
when  such  loads  are  compressive  in  nature  which  might  result  in 
local  buckling  type  failure  of  a  part  of  the  laminate  adjacent  to 
the  disbond.  Under  transverse  shear  catastrophic  and  slow  growth 
of  disbonds  may  occur  due  to  static  and  cycling  loading,  respec¬ 
tively.  Therefore,  in  this  study  attention  was  restricted  to 
growth  of  disbonds  and  failure  under  compression  and  shear  envi¬ 
ronments  . 

Development  of  new  NDT  techniques  and  modification  or  improve¬ 
ment  of  existing  ones  are  necessary  for  successful  application  of 
NDE  methodologies  to  composite  structures.  Measurements  of  stor¬ 
age  and  loss  moduli  (dynamic)  and  their  use  in  quantifying  damages 
accumulated  due  to  fatigue  and  other  environmental  factors  appear 
to  have  a  good  potential  in  that  direction.  For  this  reason  some 
exploratory  studies  were  directed  in  understanding  the  phenomenon 
of  ultrasonic  wave  propagation  in  neat  resin  and  -aminates.  Exper¬ 
imental  investigations  were  aimed  at  developing  methods  of  measure¬ 
ment  and  analysis  of  data  to  obtain  physically  meaningful  quanti¬ 
ties.  Since  matrix  material  properties  are  usually  susceptible  to 
accumulated  damages,  analytical  studies  were  aimed  at  developing 
methods  for  assessing  damages  from  measured  composite  response. 
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NDE  METHODOLOGY  TO  ASSESS  FLAW  CRITICALITY 


Nondestructive  inspection  focuses  upon  detection  of  flaws  and 
quantifying  them  in  terms  of  t^air  size  and  magnitude,  while  non¬ 
destructive  evaluation  techniques  are  needed  to  assess  their  criti¬ 
cality.  Allowable  limits  of  flaw  size  or  magnitude  must  be  deter¬ 
mined  to  establish  inspection  intervals  for  in-service  flawed 
structural  components  and  criteria  f<3£  judging  the  necessity  of 
mandatory  repair. 

For  the  interlaminar  disbonds  nondestructive  evaluation  is 
based  on  a  static/fatigue  failure  analysis.  Criticality  of  dis¬ 
bonds  under  static  transverse  shear  is  defined  in  terms  of  crack 
propagation  and  is  analyzed  by  principles  and  methods  of  fracture 
mechanics.  For  cyclic  trans verse  shear,  use  of  a  well  known  em¬ 
pirical  methodology  is  suggested  for  quantifying  damage  growth. 
Failure  of  a  disbonded  laminate  under  compressive  loading  is  de¬ 
fined  in  terms  of  buckling  and  an  elastic  stability  analysis  is 
utilized  for  assessing  criticality.  Residual  strength  or  lifetime 
can  therefore  be  predicted  by  combining  results  obtained  by  the 
use  of  techniques  described  herein. 

Analytical  methods  for  assessing  the  relationships  connecting 
speed  and  attenuation  of  ultrasonic  waves  traveling  through  the 
thickness  of  a  laminate  with  the  properties  of  its  viscoelastic 
constituents  is  presented.  These  relationships  will  be  useful  in 
assessing  the  magnitudes  of  residual  moduli  of  the  matrix  material 
from  nondestructive  wave  propagation  studies  and  hence  effects  of 
accumulated  damages,  if  any.  Experimental  NDT  methods  and  descrip¬ 
tions  of  mechanical  tests  conducted  for  verifying  the  adequacy  and 
accuracy  of  analytical  evaluation  methods  and  other  tools  are  de¬ 
scribed.  Results  of  correlation  studies  are  presented  and  dis¬ 
cussed  which  demonstrate  the  usefulness  of  the  NDE  methodology. 
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ANALYTICAL  METHODS 

PROPAGATION  OF  DISBOND  IN  STATIC  SHEAR 

Individual  layers  in  laminated  composites  commonly  used  in 
practice  are  unidirectionally  reinforced  materials  containing  fi¬ 
bers  arranged  in  a  random  array.  For  the  purpose  of  useful  stress 
analysis  such  unidirectional  composites  can  be  considered  as  ani¬ 
sotropic  materials  which  behave  elastically  until  failure.  Prin¬ 
ciples  and  mathematical  methods  employe i  in  linear  elastic  frac¬ 
ture  mechanics  are  utilized  here  for  analysis  of  stresses  near  a 
delamination  in  an  arbitrarily  laminated  structure  under  a  state 
of  plane  deformations.  In  practice  delamination  or  disbond  type 
defects  in  aircraft  components  are  usually  birth  defects  with 
two-dimensional  planar  form.  Such  defects  may  also  originate 
from  foreign  object  damage  or  various  other  reasons.  Three-di¬ 
mensional  elasticity  solutions  employing  complex  mathematical 
tools  or  finite  element  methods  must  be  obtained  for  analyzing 
stress  states  near  such  defects.  However,  2-D  defects  cannot 
propagate  in  a  self  similar  fashion  and  methods  for  studying 
their  growth  are  yet  to  be  developed,  even  for  isotropic  mate¬ 
rials.  For  design  purposes  criteria  for  one-dimensional  flaw 
growth  are  usually  employed  under  such  circumstances  which  usually 
yield  conservative  estimates  of  critical  loads.  For  a  laminated 
composite  the  state  of  the  art  is  less  advanced  because  very  few 
attempts  have  been  made  to  perform  rigorous  stress  analysis  even 
under  the  framework  of  two-dimensional  anisotropic  elasticity.  An 
approximate  strength  of  materials  approach  has  been  used  in  a  pre¬ 
vious  study  (ref.  1) .  Although  finite  element  methods  have  been 
employed  for  studying  stress  states  near  defects  in  laminated 
structures  (refs.  2,3),  it  is  known  that  such  numerical  solutions 
are  sensitive  to  various  parameters  employed  in  the  solution  (see 
ref.  4) . 

In  this  study  an  exact  solution  of  the  two-dimensional  elas¬ 
ticity  problem  of  a  simply  supported  laminated  beam  or  plate  (fig. 
1)  with  defects  is  obtained.  Solutions  to  the  governing  differen¬ 
tial  equations  for  each  layer  are  chosen  in  the  form  of  infinite 
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series,  as  outlined  in  Appendix  A-l .  Two  unknown  functions  are  in¬ 
troduced  which  characterize  the  displacement  discontinuities  at  the 
de laminations.  By  employing  a  suited) le  stiffness  formulation,  dis¬ 
placements  and  stresses  corresponding  to  each  term  of  the  series 
(harmonic  m)  are  expressed  in  terms  of  applied  load  and  the  unknown 
functions.  /  the  use  of  an  asymptotic  solution  as  m  •*>  «  and  other 
algebraic  raa  lipulation  the  problem  is  reduced  to  the  solution  of  a 
coupled  paii  of  singular  integral  equations.  The  stress  singulari¬ 
ties  at  the  tips  of  a  delamination  are  characterized.  Numerical 
solutions  are  obtained  for  thick  laminated  beams  used  in  the  test 
program.  The  parts  of  the  beam  which  contain  the  delaminations  are 
under  a  state  of  transverse  shear.  Such  a  shear  environment  alone 
is  crucial  for  aircraft  components  since  membrane  states  of  stress 
are  not  affected  by  the  presence  of  delaminations. 

Strain  energy  release  rates  for  growth  of  interlaminar  defects 
in  beam,  plate  or  shell  type  structures  have  been,  obtained  in  vari-. 
ous  studies  (refs.  5-8}  by  the  use  of  beam,  plate  or  shell  theories. 
Employing  similar  techniques  the  delaminated  beam  was  considered  in 
reference  1.  Results  of  the  present  solution  are  compared  later 
with  analytical  solutions  and  experimental  results  reported  in  ref¬ 
erence  1  and  with  test  data  obtained  in  this  program. 

As  indicated  in  Appendix  A-l,  the  coupled  pair  of  integral 
equations  are  solved  by  a  collocation  ‘method  to  determine  the  un¬ 
known  functions  ^  and  <J>2  at  a  discrete  number  of  points  over  the 
interval  [-1,1]  and  stress  intensity  factors  at  x=c  and  d  corres¬ 
ponding  to  mode  T  and  mode  II  are  evaluated  as: 

KId  =  ♦1(l)/4/2 

Klc  *  H11  ♦ 

KIId  *  H22  ^  *2<”)/4/2 

KHc  "  -H22  ^  V-1)/4/2 


(1) 
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where  H®^  depend  on  the  properties  of  the  layers  adjacent  to  the 
disbond.  In  the  test  program  the  disbonds  are  located  between  two 
similar  layers  and  in  such  cases  H^-0  and  the  stress  singularities 
at  the  tips  of  the  disbond  are  of  inverse  square  root  type.  Equa¬ 
tion  (1)  given  above  has  been  derived  based  on  this  assumption. 

When  the  disbond  is  located  between  two  dissimilar  layers  the 
singularity  is  no  longer  of  inverse  square  root  type,  but  modified 
expressions  characterizing  the  stress  distribution  near  the  tips 
of  a  delamination  can  be  derived  without  much  complication  basod 
on  the  methodology  described  in  Appendix  A-l.  However,  a  differ¬ 
ent  numerical  procedure  must  be  utilized  under  such  circumstances 
for  solving  the  integral  equation.  For  inverse  square  root  type 
singularity  the  strain  energy  release  rate  is  given  by: 


(2) 


Catastrophic  propagation  of  the  disbond  will  occur  when  G  reaches 

the  critical  value  G  ,  which  can  be  considered  as  a  material  con- 

c 

stant.  The  results  of  the  experimental  program  are  utilized  to 
evaluate  Gc.  In  a  shear  environment  Kj  is  extremely  small  com¬ 
pared  to  KIX  and  catastrophic  failure  may  also  be  characterized  by 
a  fracture  toughness. 


GROWTH  OF  DISBOND  IN  CYCLIC  SHEAR  ENVIRONMENT 

Fatigue  crack  propagation  is  a  complex  phenomenon  which  de¬ 
pends  on  the  nature  of  applied  loads,  fatigue  and  mechanical  pro¬ 
perties  of  the  material,  micros true tural  process  zones  in  the 
immediate  vicinity  of  the  crack  tip  and  various  other  factors. 
Other  phenomena  which  may  have  significant  influence  on  growth  of 
a  flaw  are  crack  blunting,  crack  closure,  residual  stresses,  etc. 
A  realistic  analytical  model  for  crack  propagation  has  yet  to  be 
developed.  Some  studies  in  this  direction  (refs.  9-12)  have  been 
attempted  for  ductile  metals  with  limited  success.  Since  micro- 
structural  processes  near  a  disbond  in  composites  are  not  yet  as 
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clear  as  these  in  metals  and  not  much  is  known  about  fatigue  pro¬ 
perties  of  complex,  anisotropic  composite  laminae  and  laminates,  a 
practical  approach  to  study  propagation  of  disbonds  in  laminated 
composites  appears  to  be  the  use  of  the  well  known  semi-empirical 
crack  growth  law,  i.e.: 


if  -  cii«n 


•  (3) 


where 

a  ■  disbond  or  crack  length 
N  =  number  of  cycles 

AK  *  the  stress  intensity  factor  range  ■  r\  /a  (Stti_v.-Stti.  _) 

Smnv,Sm.  ■  the  maximum  and  minimum  value  of  applied  load 
max  min 

or  stress 

n  s  a  factor  depending  on  geometry  and  other  variables  in¬ 
fluencing  stress  intensity  factor 
C,n  »  empirical  constant  and  exponent  obtained  from  curve  fit 
to  experimental  data. 

In  this  study,  K^,  the  mode  II  stress  intensity  factor  is 
utilized  as  K  and  a  is  calculated  based  procedures  outlined  in 
the  preceding  subsection  and  Appendix  A-l,  which  yield  stress  in¬ 
tensity  factors  based  on  the  exact  solution  of  the  two-dimensional 
anisotropic  elasticity  problem  of  the  laminated  beam  considering 
the  stacking  sequence  of  the  laminate. 

BUCKLING  FAILURE  OF  COMPRESSION  SKIN  CONTAINING 
A  DELAMINATION 

When  a  laminate  containing  a  disbond  is  subjected  to  compres¬ 
sive  loading,  buckling  type  failure  may  result  in  a  significant 
reduction  in  strength.  Such  failures  are  likely  to  be  most  criti¬ 
cal  in  the  compression  skin  of  a  sandwich  beam  or  plate.  In  this 
study  attention  is  restricted  to  the  one-dimensional  beam  type 
structure  containing  a  single  disbond.  For  the  purpose  of  analysis 
the  compression  skin  is  considered  to  be  supported  on  an  elastic 
foundation  providing  constraints  in  shear  as  well  as  extensional 
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loadings.  The  skin  is  modeled  as  an  assemblage  of  four  beam  ele¬ 
ments  (see  fig.  2) ,  namely: 

1.  Two  beams  (1,2)  of  finite  length  equal  to  the  length  of 
the  disbond,  one  above  and  the  other  below  it;  and 

2.  Two  semi-infinite  beams  (3,4)  connected  to  the  finite 
beams  at  the  tips  of  the  disbond. 


Differential  equations  governing  the  deformations  of  the  com¬ 
pressed  laminated  beam  elements  resting  on  elastic  foundation  in¬ 
clude  the  effects  of  shear  deformation.  These  equations  are 
solved  and  stiffness  matrices  for  the  elements  are  derived  for 
symmetric  deformation  patterns,  since  the  lowest  value  of  criti¬ 
cal  load  is  desired.  Determinants  of  the  global  stiffness  matrix 
are  calculated  for  a  set  of  prescribed  compressive  loads  and  the 
critical  load  is  determined  from  the  condition  that  the  determin¬ 
ant  is  equal  to  zero.  A  trial  and  error  procedure  is  utilized  to 
determine  the  eigenvalue  up  to  the  desired  accuracy.  Details  of 
the  formulation  are  jiven  in  Appendix  A-2.  Sandwich  beam  speci¬ 
mens  (fig.  3)  with  varying  disbond  lengths  were  designed  to  test 
the  validity  of  the  flaw  criticality  criterion  determined 
from  analysis.  Details  of  the  test  specimens  can  be  found  in  the 
following  section.  Maximum  disbond  length  was  £«38.1  mm  in  a  beam 
of  span  L-508  mm.  Therefore,  the  assumption  that  beams  3  and  4 
are  semi- infinite  in  length  appears  to  be  reasonable.  The  thick¬ 
ness  of  the  12  ply  compression  skin  (Q/±45/*45/0) _  with  the  dis- 
bonds  located  in  the  mid-plane  was  small  (of  the  order  of  1.5  mm) 
as  compared  to  the  depth  of  the  aluminum  honeycomb  core  (Hc»38.1 
mm).  Therefore,  the  skin  can  be  considered  to  be  under  a  state 
of  uniform  compressive  strain.  Compressive  forces  on  the  four 
beam  elements,  however,  differ  in  magnitude.  Axial  and  bending 
stiffnesses  of  the  beam  elements  are  computed  using  laminate  anal¬ 
ysis,  Shear  stiffnesses  are  calculated  on  the  assumption  of  con¬ 
stant  shear  stress  in  the  layers,  i.e.,  l/Cg5«E  vVc^,.,  V*  being 


'55' 


the  volume  fraction  of  layer  i.  In  addition,  a  shear  correction 
factor  k55  equal  to  5/6  is  used  to  obtain  the  effective  shear 
stiffness.  For  the  purpose  of  computation  of  foundation  modulus 
the  core  is  assumed  to  be  fixed  at  the  tension  flange. 
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TRANSVERSE  WAVE  PROPAGATION  AND  COMPLEX  MODULI 

An  important  nondestructive  evaluation  test  consists  of 
sending  a  pulse  transversely  through  the  laminate.  Compari&on  of 
transmission/reflection  data  for  unflawed  and  flawed  laminates 
may  then  serve  to  uncover  characteristic  flaws  in  this  manner. 

The  following  discussion  will  be  concerned  with  simple  transverse 
wave  propagation  through  a  laminate  and  of  the  complex  moduli 
associated  with  matrix  deunping  and  the  waves  considered.  Propa¬ 
gation  of  elastic  waves  in  a  laminate  is  duscussed  in  Appendix 
A- 3. 

A  polymer  fiber  composite  will  exhibit  viscoelastic  effects 
if  the  polymeric  matrix  viscoelasticity  is  significant.  If 
sinusoidal  waves  are  transmitted  transversely  to  the  laminate 
then  it  follows  from  (A- 3. 10)  and  standard  viscoelastic  wave 
theory  (see  ref.  13)  that  the  wave  speed  is: 


and  the  amplitude  attenuation  is: 


where  &  and  &Tare  the  complex  counterparts  of  the  elastic  constant:; 
k  and  G^. 

$:  (ico)  *  k'  (<o)  +  ik"  (w)  (6) 

&  ( i<u )  »  G  (to)  +  iG  («) 

Here  <o  is  the  frequency,  prime  and  double  prime  denote  real  and 
imaginary  part,  respectively,  and  i  ■ 
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The  experimental  determination  of  the  complex  moduli  k  and 
entering  into  (4,5)  is  a  matter  of  considerable  difficulty. 
Fortunately,  however,  these  complex  moduli  can  be  computed  on  the 
baais  of  matrix  viscoelastic  and  fiber  elastic  properties  by 
methods  given  in  reference  14.  Accordingly  the  complex  effective 
moduli  and  are  given  by  the  expressions: 


Kikt*K)vn  *  kf(Va,n>vl 
<lt£+an,WB  +  (Wvf 


(7) 


/V 


G  + 
m 


+  2^m 
m  m 


G-f-&  2&(Jc+<S)  m 

Tt  m  m  m  m 


(8) 


Here: 

m  -  Index  indicating  matrix 
f  -  Index  indicating  fiber 
v  -  Volume  fraction. 


The  isotropic  matrix  complex  plane  strain  bulk  modulus  can 


be  written: 

k  »  K  + 

i 

i' 

(a) 

m  m 

3" 

m 

k '  "  Km  + 

tt  m 

1 

3 

Gm 

(b) 

n 

k  * 

m 

1 

T 

ll 

Gu 

(c) 

where  K  is  the  three-dimensional  elastic  bulk  modulus  and  the 
in 

usual  assumption  has  been  made  that  matrix  dilatational  visco¬ 
elasticity  can  be  neglected. 
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The  matrix  complex  shear  modulus  is  written  in  the  form 

(6)  : 


&(i<u)  -  g' (<o)  +  1g"  (co  >  (a) 

m  m  m 

(10) 

tan  »  g"/g'  (b) 


where  the  last  expression  is  the  loss  tangent. 

1.  is  now  necessary  to  find  the  real  and  imaginary  parts  of 
(8).  This  tedious  undertaking  is  greatly  facilitated  by  exploi¬ 
tation  of  the  fact  that  (10)  is  generally  a  small  number  (of 
order  .05).  It  can  be  shown  (ref.  15)  that  in  that  case: 


i  i 

k  *  k(k.G) 
m  m 


k"  -  k" 

3Gm 


(11) 


&(k’  G‘) 
m  m 


3G 


(12) 


m 


3G, 


m 


It  follows  that: 


,  k'(k_+G_)v  +  k.(k  +g’)v. 

(k,+Gjv  +  (k ' +G ' )  V- 
r  m  m  m  m  £ 


(a) 
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This  implies  chat  a  longitudinal  v.a  vs  propagating  transversely 
through  the  laminate  has  velocity  given  by  (15)  and  its  amplitude 
after  propagating  distance  d  is  attenuated  by  the  amount  exp  (>oad) . 

The  wave  speed  (15)  and  the  attenuation  (16)  can  be  deter¬ 
mined  experimentally  by  propagating  a  longitudinal  wave  normal  to 
the  laminate.  Examination  of  (13),  (14)  which  enter  into  (15), 

9  H 

(16)  shows  that  the  real  and  imaginary  parts  G_  and  G_  off  the  ma¬ 
trix  shear  modulus  can  be  determined  from  (15),  (16)  if  all  other 
quantities  entering  are  known.  These  are  fiber  elastic  proper¬ 
ties  G__  and  k«,  fiber  and  matrix  volume  fractions  v.  and  v_  and 

•  I  X  ^  X  Itl 

the  real  part  km.  Recalling  (9b) ,  it  is  seen  that  for  the  last 

quantity  it  is  merely  necessary  to  know  the  frequency  independent 

three-dimensional  matrix  bulk  modulus  K  .  It  is  thus  seen  that 

m 

normal  wave  propagation  can  serve  to  determine  in  situ  values  of 

I  II 

G  and  G  . 
m  m 
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EXPERIMENTAL  PROGRAM 
SHEAR  TESTS  ON  THICK  DISBONDED  LAMINATES 

The  thick  beam  laminates  were  tested  in  three  point  bend¬ 
ing  with  a  fixture  specially  manufactured  for  these  te3ts. 

The  specimen  dimensions  were  254  mmx  25.4  mm  x  8.9  mm  (10  in.  x 
1  in.  x  .35  in.).  The  support  span  was  152.4  mm  (6  in.).  The 
center  support  consisted  of  two  11.13  mm  (0.438  in.)  diameter 
cylinders  which  were  clamped  to  each  side  of  the  specimen.  This 
support  was  attached  to  wedge  aiction  friction  grips  which  were 
connected  to  the  load  cell.  The  end  supports  were  also  11.13  mm 
(0.438  in.)  cylinders  which  only  contacted  one  aide  of  the  specimen 
(fig.  4). 

The  static  tests  were  performed  on  an  Instron  (Model  TTC) 
static  testing  machine.  The  specimens  were  loaded  at  a  cross 
head  rate  of  .05  cm/min.  The  failure  load  was  recorded  from 
the  chart  recorder  as  the  maximum  load  experienced  prior  to  fail¬ 
ure.  The  fatigue  tests  were  performed  in  an  Instron  Servo-hy¬ 
draulic  fatigue  testing  machine.  The  load  was  varied  in  a  sinusoid 
ut  a  frequency  of  10  HZ  with  the  minimum  load  equal  to  ten  percent 
of  the  maximum  load.  Several  maximum  load  levels  were  tested. 

Two  different  graphite-epoxy  laminates  were  tested.  The 
[  (0/90)g012) a  laminate  had  two  types  of  implanted  defects  (fig.  5). 

The  25.4  mm  (1.0  in.)  square  defect  was  oriented  such  that  the 
delamination  was  perpendicular  to  the  longitudinal  axis  of  the 
specimen.  The  25.4  mm  (1.0  in.)  diamond  defect  was  oriented  at 
45°  to  the  specimen  axis  so  that  the  delamination  came  to  a  point. 
Each  specimen  contained  two  implanted  defects  located  midway  be¬ 
tween  the  center  and  end  supports.  These  specimens  were  pre¬ 
cracked  by  clamping  the  specimens  through  the  thickness  to  arrest 
crack  growth  and  then  loaded  until  cracking  was  detected  by 
acoustic  emissions.  The  specimens  were  then  ultrasonically  C- 
scanned  to  determine  the  precrack  length.  These  lengths  are 
given  in  Table  1.  The  specimens  were  then  statically  tested 
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without  clamps.  The  failure  loads  for  these  tests  are  also  given 
in  Table  1. 

The  t (0^/+452/+452/0^ ) g] g  laminate  was  implanted  with  diamond 

defects  only.  The  defect  lengths  were  19.05  mm  (0.75  in.),  25.4  mm 
(1.0  in.),  and  31.75  mm  (1.25  in.).  These  samples  were  tested  statically 
with  no  precracking.  The  static  test  results  are  given  in  Table  2.  The 
fatigue  tests  were  conducted  at  several  different  S-levels. 

Samples  with  0.75  in.  and  1.0  in.  defects  were  fatigued  at  S- 
levels  of  0.4,  0.5  and  0.6.  The  samples  with  1.25  in.  defects 
were  fatigued  at  S-levels  of  0.3,  0.4,  0.5  and  0.75.  The  number 
of  cycles  to  failure  are  given  in  Table  3.  Typical  flaw  prop¬ 
agation  results  are  presented  in  fig.  6-17. 

COMPRESSIVE  BUCKLING  OF  DISBONJED  LAMINATES 

The  compressive  buckling  properties  of  an  AS-3 501-6  graphite- 
epoxy  laminate  were  determined  by  four  point  bending  of  a  sand¬ 
wich  beam  specimen.  The  specimens  were  constructed  from  24  lb. 

Hexel  aluminum  honeycomb  adhesively  bonded  to  two  12  ply  graphite- 
epoxy  faces.  The  honeycomb  thickness  was  38.1  mm  (1.5  in.)  and 
the  sample  width  was  25.4  mm  (1  in.).  The  overall  specimen  length 
was  559  mm  (22  in.).  Figure  18  gives  the  specimen  and  loading 
geometry. 

The  laminate  which  was  tested  was  the  [0/+45-/0]  .  Two  layers 

s 

of  1  mil  teflon  were  embedded  in  the  laminate  to  simulate  a  de¬ 
lamination.  The  near  surface  defect  was  implanted  between  the 
third  and  fourth  plies  while  the  center  disbond  was  implanted 
between  the  two  middle  plies.  The  length  of  these  rectangular 
disbonds  was  varied  from  12.7  mm  (0.5  in.)  to  38.1  mm  (1.5  in.). 
Specimens  with  no  implanted  defects  were  also  tested. 

For  each  tests  the  longitudinal  strain  on  the  compressive 
face  (ec)  was  plotted  versus  the  strain  on  the  tensile  face  (eT) 

on  an  X-Y  recorder.  Since  the  strain  gages  were  located  directly 
over  the  implanted  defect,  it  can  be  seen  that  when  the  specimen 
buckled  at  the  defect,  the  compressive  strain  was  reduced.  As 
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the  load  is  increased,  the  buckled  region  "bows  out"  and  the 
compressive  strain  is  reduced  accordingly.  For  the  samples  with 
19.05  mm  (0.75  in.)  debonds  and  the  samples  with  no  defects,  the 
compressive  strain  was  also  plotted  against  the  load  which  was 
applied  to  the  beam.  The  shape  of  these  curves  is  identical  to 
the  strain  versus  strain  curves.  The  buckling  load  was  taken  to 
be  the  load  at  which  the  maximum  compressive  strain  was  attained. 
This  load  was  taken  from  the  Instron  chart  for  the  samples  where 
load  versus  compressive  strain  was  not  plotted.  The  stress  at 
buckling  was  calculated  from  beam  theory  to  be: 


where  P  is  the  applied  load  at  buckling,  t  is  the  laminate  thick¬ 
ness  (1.524  mm  «  .06  in.),  and  w  is  the  specimen  width  (25.4  mm  * 
1.0  in.).  The  experimental  results  are  summarized  in  Table  4. 

The  results  for  the  specimens  with  no  defects  show  that 
both  the  ec  versus  eT  and  ec  versus  load  curves  are  linear  to 
failure  (fig.  19). 

The  behavior  of  the  samples  with  defects  in  the  center  of 
the  laminate  was  highly  dependent  upon  the  length  of  the  im¬ 
planted  defect.  For  the  smaller  size  defects  (0.5  in.  and  some 
of  the  0.75  in.  samples),  the  curves  were  linear  to  failure 
(fig,  20,21).  The  failure  stress  of  these  samples  was  consider¬ 
ably  less  than  that  for  the  samples  with  no  defects.  This  means 
that  the  presence  of  the  smaller  defects  does  affect  strength 
but  this  behavior  can  not  be  detected  from  the  recorded  strains. 
For  the  samples  with  larger  defects  (some  of  the  0.75  in.  and 
all  1.0  in.  samples),  the  e_  Versus  and  e„  versus  load  curves 
were  linear  until  the  samples  begem  to  buckle  at  the  debond 
(figs.  22,23).  When  this  occurred,  the  compressive  strain  was 
reduced  as  the  area  over  the  debond  "bows  out."  For  the  samples 
with  the  1.5  in.  defects  in  the  center  of  the  laminate,  the  be¬ 
havior  was  the  same  as  the  1.0  in.  defrtts  except  for  an  initial 
non-linearity  (fig.  24).  This  is  due  to  the  fact  that  for  these 
samples,  the  tension  side  of  the  beam  was  prebuckled  due  to  re¬ 
sidual  thermal  stresses  in  the  1.5  in.  near  surface  defects. 

Thus  there  is  an  initial  reduction  in  tensile  strain  until  this 
prebuckle  "flattens  out." 
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The  samples  with  the  near  surface  defects  were  prebuckled 
due  to  thermal  residual  stresses  from  the  curing  or  bonding 
operations.  As  a  result,  the  application  of  a  compressive  load 
to  these  samples  causes  the  buckle  to  "bow  out"  immediately. 

Thus  the  compressive  strain  becomes  negative  (or  becomes  a  ten¬ 
sile  strain).  The  1.0  in.  near  surface  defect  samples  showed  a 
small  increase  in  compressive  strain  before  buckling  occurred 
(fig.  25).  The  stress  at  buckling  was  very  small.  The  1.5  in. 
defect  samples  were  buckled  prior  to  testing  and  as  a  result  the 
comprssive  strain  decreased  linearly  with  the  tensile  strain  of 
the  tensile  aide  of  the  beam  (fig.  26).  Therefore,  no  buckling 
stress  was  measured  for  these  samples. 

In  summary,  it  can  be  seen  for  the  center  defects,  the 
compressive  stress  at  failure  decreases  with  defect  length 
(fig.  27).  The  stress-strain  response  for  smaller  defects  is 
linear  until  failure.  For  the  larger  defects,  the  material  over 
the  defect  buckles,  causing  a  reduction  in  conqpressive  strain. 

The  prebuckling  phenomenon  of  the  near  surface' defect  samples  caused 
the  samples  to  buckle  immediately  after  the  application  of  load. 

For  all  samples  with  implanted  defects,  failure  was  due  to  buck¬ 
ling  of  the  material  above  the  defect  which  resulted  in  a  de- 
lamination  extending  from  the  implanted  defects  (figs.  28,  29,  30). 

MECHANICAL  CHARACTERIZATION  FOR  MOISTURE  CONDITIONING  EFFECTS 

In  an  effort  to  develop  environmental  degradation  in  composite 
laminates  for  subsequent  nondestructive  evaluation,  test  sample 
laminates  of  [+45] 2g  configuration  were  exposed  to  100  percent 
relative  humidity  conditions  (vapors  of  boiling  water)  for  a 
period  of  four  months.  The  samples  ware  removed  and  tested  in 
a  dry  condition  in  order  to  establish  mechanical  property  de¬ 
gradation.  Table  5  lists  the  results  of  these  tests. 

Samples  of  the  neat  resin,  3501-6  were  also  subjected  to  the 
environmental  exposure  discussed  above.  No  mechanical  characteriza¬ 
tion  of  the  neat  resin  was  undertaken. 
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ULTRASONIC  WAV^  PROPAGATION  STUDIES 

Variable  Incidence  Ultrasonic  Test  Facility 

To  evaluate  the  ultrasonic  material  elastic  constraints 
other  than  C33,  it  is  necessary  tc  measure  both  the  longitudinal 
and  shear  velocities  of  sound  in  the  material.  Shear  waves  are 
produced  in  a  sample  when  it  is  inclined  at  an  angle  not  normal 
to  the  incoming  sound.  The  mode  conversion  resulting  from  non¬ 
normal  incidence  produces  waves  with  both  shear  and  longitudinal  po¬ 
larization  as  shown  in  fig.  31.  Due  to  the  translation  of  the 
shear  and  longitudinal  velocity  components ,  it  is  necessary  to 
be  capable  of  translating  either  -die  transmitting  or  receiving 
transducer  90°  to  the  transmitted  sound  as  well  as  being  capable 
of  rotating  both  the  specimen  and  the  receiving  transducer. 

In  order  to  perform  through  ‘transmission  ultrasonic  measure¬ 
ments  in  both  the  normal  incidence  and  non-normal  incidence 
mode,  a  test  tank  was  constructed.  The  tank  consists  of  two 
concentric  turntables  with  a  motor  drive  and  position  encoder 
for  each.  The  center  turntable  has  a  vise  for  holding  the 
sample  and  the  outer  turntable  has  a  mount  for  the  receiving 
transducer.  The  transmitting  transducer  is  mounted  on  the  tank 
base.  The  tank  is  constructed  of  anodized  aluminum  and  plexi- 
glas..  The  motor  drive  units  are  made  of  steel  and  are  attached 
to  the  turntables  via  a  worm  drive  system.  The  encoders  are 
attached  to  the  turntables  through  a  cable  mechanism  which  drives 
a  screw  attached  to  two  15  turn  precission  potentiometers.  The 
turntable  position  is  given  by  a  voltage  output  which  ranges 
from  0.00  to  3.60  volts. 

The  transmitting  and  receiving  transducers  *re  connected 
to  a  panametrics  ultrasonic  analyzer  and  the  output  is  monitored 
with  a  Nicolet  digital  oscilloscope  set  for  a  sampling  rate  of 
20  MHz.  The  20  MHz  sampling  rate  implies  that  only  signals  up 
to  10  MHz  may  be  observed  without  aliasing; therefore,  in  order 
to  satisfy  the  Nyguist  criteria,  a  low  pass  filter  with  a  10  MHz 
cut-off  frequency  is  inserted  between  the  output  of  the  ultra¬ 
sonic  analyzer  and  the  oscilloscope  input.  The  oscilloscope  is 
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capable  of  outputting  information  in  digital  form;  however, 
before  this  information  could  be  analyzed,  an  interface  had  to 
be  built  to  connect  the  digital  oscilloscope  to  an  LSI  11  Digital 
computer.  The  interface  uses  standard  TTL  components  to  form  a 
multiplexed  seven  channel  input  output  port.  Each  channel  has 
sixteen  input  and  sixteen  output  bits.  The  output  channels  have 
both  latching  and  tristate  capability.  The  data  from  the  oscillo¬ 
scope  is  transmitted  through  the  LSI  11  computer  to  a  Burrows  B7700 
computer  over  a  telephone  line  for  analysis  and  storage.  After 
the  data  is  stored,  it  can  be  analyzed  using  FFT  techniques  an' 
a  frequency  spectrum  of  the  ultrasonic  waveform  can  be  produced 
In  addition  to  analyzing  the  total  waveform,  an;-  part  of  the  data 
may  be  analyzed  at  any  time  without  the  need  of  performing  the 
test  again. 

Software  Development 

LSI  11  machine  language  programs  had  to  be  developed  in 
order  to  control  the  interface  between  the  digital  oscillascope 
and  the  computer.  This  software  has  the  capability  of  initial¬ 
izing  the  interface  and  transmitting  any  portion  of  the  ultra¬ 
sonic  waveform  through  a  modem  connected  to  a  telephone  line. 

The  program  fills  a  file  with  data  after  the  file  is  initially 
created  by  the  operator.  The  operator  can  communicate  with  the 
computer  through  a  Digital  LA36  line  printing  terminal  or  through 
a  Tektronix  4006-1  graphics  terminal.  The  graphics  terminal  is 
used  to  plot  results  of  the  FFT  routines.  This  data  transmission 
program  is  called  NICDATA 

In  order  to  determine  the  frequency  spectrum  of  the  wave¬ 
form,  a  program  called  NICTRUM  was  written  in  fortran  wnich  runs 
on  the  Burrows  B7700  computer.  The  program  takes  the  FFT  of  the 
entire  data  string  and  displays  both  the  ultrasonic  waveform  and 
the  spectrum  on  either  a  Tektronix  terminal  or  flatbed  plotter. 

In  order  to  analyze  selected  portions  of  the  total  waveform, 
a  program  called  PARTITION  was  written.  This  program  allows 
the  data  string  to  be  plotted  and  then  portions  of  the  waveform 
can  be  partitioned  for  analysis  separate  from  the  total  waveform. 
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A  program  which  takes  the  FFT  of  the  first  partition  was  written 
and  is  called  SELECTRUM  (fig.  32).  Another  program  was  written 
to  compare  spectrums  of  different  parts  of  the  waveform. 
COMPARETRUM  (fig.  33)  takes  the  FFT  of  the  first,  second  and 
third  partitions  as  well  as  the  total  waveform  and  plots  them 
on  one  graph  as  either  dotted  solid  or  dashed  line. 

PHASE  is  a  program  which  evaluates  the  phase  of  the  first 
and  second  partitioned  portions  of  the  waveform.  These  par¬ 
titioned  sections  are  normally  the  first  and  second  echos. 

Such  phase  measurements  are  necessary  for  evaluating  velocity 
of  sound  measurements  (fig.  34) . 

VELOCITY  is  a  program  which  uses  the  phase  information  to 
determine  the  velocity  of  sound  as  a  function  of  frequency  by 
the  equation  (fig.  35) 

,  .  2  do 

V  W  <j»2  (aj)H> i  (w)  (18) 

where:  d  is  the  sample  thickness 
to  is  the  angular  velocity 
<|>2(o)  is  the  phase  of  the  second  echo 
<pi(to)  is  the  phase  of  the  first  echo 


ATTENUATION  (fig.  36)  is  a  program  which  uses  information 
from  the  FFT  of  three  partitioned  portions  of  the  waveform.  The 
partitioned  portions  are  the  front  and  back  surface  echos  and 
the  second  multiple  echo.  Using  this  information  the  acoustic 
attenuation  a  can  be  evaluated  as  a  function  of  frequency  by  the 
equation 
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P 


(19) 


R^Cto),  R,  (to)  and  R^  (to)  are  the  magnitudes  of  the  Fourier  trans¬ 
forms  of  the  three  partitioned  waveforms. 

A  program  was  also  written  to  list  the  data  stored  from 
the  oscilloscope  in  a  tabular  form  along  with  thickness  and  par¬ 
tition  information.  This  program  was  called  NICTABLE. 


-22- 


NADC-77278-30 


Normal  Incidence  Ultrasonic  Pulse-Echo  Measurements 

Normal  incidence  pulse-echo  measurements  were  made  on  +45° 
graphite  samples,  fifty  percent  of  which  had  been  moisture 
conditioned  in  a  90°  C  steam  bath  until  they  were  saturated 
with  moisture.  Normal  incidence  pulse-echo  measurements  were 
also  made  on  non-moisture  conditioned  neat  resin  samples.  The 
resin  type  was  that  used  in  the  AS3501-6  graphite/ epoxy  system. 
The  acoustic  attenuation  coefficient  a  was  evaluated  by  measur¬ 
ing  a  single  sample  thickness  and  the  amplitudes  of  the  front 
surface,  back  surface  and  second  multiple  echos.  The  acoustic 
attenuation  is  found  by  using  the  equation 

1  ^2 

°  =  "  ln  (Rj  "  (20) 


where : 


is  the  linear  amplitude  attenuation  coefficient  in  mm" 

is  the  sample  thickness  in  mm 

is  the  amplitude  of  the  front  surface  echo 

is  the  amplitude  of  the  back  surface  echo 

is  the  amplitude  of  the  second  multiple  echo 


The  results  of  the  pulse  echo  attenuation  measurements  are 
summarized  in  Table  6.  The  high  degree  of  scatter  in  the  dry 
+45°  graphite  samples  is  due  to  the  thin  sample  thickness  and 
the  difficulties  associated  with  separating  echos  in  thin 
laminates.  Additional  measurements  were  made  on  thick  beam 
laminates  and  the  resulting  average  attenuation  for  dry  graphite 
epoxy  was  evaluated  to  be  .0125  mm-*.  The  average  attenuation 
in  the  moisture  saturated  graphite  epoxy  was  evaluated  as  .0543 
and  the  neat  resin  had  an  average  acoustic  attenuation  of  .113  mm-1. 
Average  velocities  of  sound  were  also  measured  for  normal  incidence 
sound  on  graphite  epoxy  and  neat  resin.  The  values  are  3190  m/sec 
for  graphite  and  2750  m/sec  for  the  resin  used  in  the  AS3501-6 
graphite  epoxy  system. 
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ANALYTICAL/EXPERIMENTAL  DATA  CORRELATION 

STATIC  SHEAR  OF  THICK  LAMINATED  BEAMS 

Numerical  results  for  a  t  (04/+452/+452/04)  ]  laminate 

with  64  plies  and  H/&  ratio  of  .05872  considered  in  reference  1 
obtained  by  the  use  of  the  methods  presented.  Two  solu— 
tions  obtained  are  discussed  below. 

1.  Laminate  stacking  sequence  is  retained  but  (+452/+452) 
layers  are  chosen  as  one  material  with  average  pro¬ 
perties.  Therefore  the  beam  has  nine  layers  of  ap¬ 
propriate  thicknesses  with  delamination  at  the  mid¬ 
plane  of  the  0°  layer  at  the  center. 

2.  Since  average  properties  (or  equivalent  elastic 
moduli;  are  often  used  in  mechanics  of  composite 
materials  a  second  solution  was  obtained  by  using 
average  properties  for  the  whole  structure,  i.e.  a 
single  layer  with  disbond  at  the  mid-plane. 

All  elastic  properties  are  listed  below  and  are  nondimen- 
sionalized  by  a  factor  EQ=6.895  GPa 


Layer* 

C11/E0 

C33/£0 

C13/E0 

C55/Ei 

1)  c 

18.24 

1.46 

.409 

.841 

2)  (+  2/+452) 

2.88 

1.50 

.134 

.687 

3)  Average 

10.92 

1.48 

.2 

.756 

These  va.  ies  are  representative  of  those  measured  from  exper¬ 
iments  (reference  1).  For  averaging  the  properties,  constant 
strain  .  J  constant  stress  hypotheses  were  used  for  computing 
stiffnesses  parallel  and  perpendicular  to  laminations.  39 
terms  of  the  series  and  20  collocation  points  were  found  to 
yield  convergent  solutions  for  all  practical  purposes. 

Nondimensionalized  stress  intensity  factors  and  strain  energy 
release  rates  from  the  two  elasticity  solutions  are  shown  in 
Table  7  for  various  values  of  a/£.  It  is  noted  that  if  the 
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beam  is  considered  as  homogeneous,  then  the  problem  is  equivalent 
to  that  with  delaminations  with  applied  shear  stress  equal  to 
1.5  P/2H  on  the  surfaces,  provided  the  delaminations  are  away 
from  the  load  and  the  supports.  For  a/1  and  a/h  approaching 
zero,  the  limiting  value  of  KIX  should  be  1.0607  P/S/2H.  Although 

KJX  is  higher  than  this  value  for  large  a/1,  it  is  smaller  for 
a/1  «  .004.  Therefore,  as  a/1  is  reduced  further,  Klx  should 

approach  the  limit  from  below.  This  trend  is  similar  to  the 
trends  reported  in  Reference  16.  Values  of  Kxx  at  both  tips 

differ  only  in  the  third  significant  figure  (given  here  i  to  two 
decimal  places).  Mode  1  intensity  factors  are  extremely  small 
and  will  not  have  any  influence  on  propagation  of  the  disbonds. 

This  is  expected  since  H/l  is  only  .05872.  Therefore,  the  con¬ 
jecture  made  in  Reference  1  regarding  effects  of  concentrated 
loading  on  retarding  crack  growth  possibility  is  incorrect. 

For  larger  disbond  lengths  such  effects  may  have  some  influence. 

One  interesting  fact  is  that  the  nondimensionalized  value 
of  KIX  from  both  theories  are  closer  for  smaller  crack  lengths 

but  the  strain  energy  release  rates  are  far  apart.  It  is  clear 
that  stresses  will  have  similar  distributions  near  the  tips  of 
the  disbonds  for  small  a /I .  The  difference  in  G  occurs  due  to 
the  difference  of  in  the  two  theories.  Percentage  difference 
in  G,  however,  decreases  as  a/1  is  increased.  Therefore  averaged 
properties  may  be  used  with  proper  caution  and  correction  to  ob¬ 
tain  reasonably  accurate  results.  Critical  loads  and  strain 
energy  rates  are  compared  with  beam  theory  solutions  and  exper¬ 
imental  data  in  Table  8  and  fig.  37.  Fig.  37  shows  that  beam 
theory  yields critical  loads,  which  are  larger  than  elasticity 
solutions.  The  difference  is  large  for  small  a/1  ratios  but 
is  gradually  reduced  as  a/1  increases.  To  some  extent  this  re¬ 
duces  the  discrepancy  in  measured  G_  values  for  two  disbond  lengths 
The  difference  in  G_  computed  from  experimental  data  by  the  use 

v 

of  beam  theory  solution  is  44%,  whereas  if  elasticity  solution 
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is  utilized,  the  difference  is  reduced  to  20%,  which  is  less  than 
the  experiments 1  scatter  for  each  disbond  length* 

Critical  values  of  strain  energy  release  rates  obtained 
by  the  use  of  results  from  tests  on  precracked  [  (0/90) ^0 s 

laminates  and  exact  elasticity  solutions  considering  stacking 
sequence  are  given  in  Table  9.  After  comparing  these  results 
with  those  in  Table  8  (exact  elasticity)  it  appears  that  values 
of  Gc  obtained  from  (i)  precracked  ((0/90)  6°12>s  specimens 

and  (ii)  l (04/+452/+452/04 ) g]  laminated  beams  without  any  pre¬ 
cracking  of  the  disbonds,  are  of  the  same  order  and  show  same 
amount  of  scatter.  Precracking,  therefore,  does  not  have  any 
significant  influence  on  catastropic  disbond  growth  under  static 
shear. 

PROPAGATION  DISBOND  IN  FATIGUE 

The  fatigue  test  samples  with  diamond  shaped  defects  could 
not  be  precracked  to  obtain  a  straight  line  crack  front.  Fatigue 
tests  were  therefore  conducted  with  the  implanted  shapes.  Lengths 
of  extended  disbonds  were  measured  from  "C"  scans  with  an  ac¬ 
curacy  of  +.6  mm  after  specified  number  of  cycles, N, at  six  loca¬ 
tions,  three  in  each  implanted  disbond  in  the  thick  beam  specimens 
as  shown  in  figure  38.  it  is  clear  from  figures  6-17  that  at 
locations  3  and  6,  both  at  mid-width  of  the  beams,  little  or  no 
extension  occurred.  The  reason  for  this  phenomenon  is  not  quite 
clear  but  may  be  attributed  to  the  fact  that  the  stress  distri¬ 
bution  near  these  corners  are  different  in  nature  as  compared  to 
that  near  the  edges  (see  reference  17).  Extension  of  disbonds 
at  locations  1,  2,  4  and  5  were  also  not  uniform,  i.e.  signif¬ 
icantly  more  extensions  occurred  at  some  locations  as  compared 
to  others.  This  behavior  is  expected  because  of  local  non-uniform¬ 
ities  in  materials  and  shapes  of  the  disbond  tips. 

Disbond  lengths  at  locations  1,  2,  4  and  5  were  plotted 
against  N  and  reasonably  smooth  curves  were  obtained.  Some 
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sample  plots  are  shown  in  figures  39-42.  Figures  40-42  show 
that  in  some  locations  a  plateau  exists  in  a  vs.  N  plots.  This 
type  of  behavior  is  usually  observed  in  metals  when  sudden  over¬ 
loads  occur  at  a  certain  point  in  a  loading  cycle.  No  apparent 
reason  was  found  for  this  behavior  from  the  testing  program. 

It  is  possible  that  material  nonhomogeneities  encountered  during 
crack  growth  (existence  of  weak  locations)  may  also  result  in 
such  behavior.  Nonhomogeneities  in  composites  are  usually  much 
more  severe  them  those  in  metals.  It  should  be  pointed  out, 
however,  that  this  type  plateau  were  observed  in  one  or  two  lo¬ 
cations  a  limited  number  of  samples,  namely  .75-5,  1.00-4,  1.00-6 
and  1.25-7. 

Slopes  Aa/AN  of  a-N  curves  were  then  calculated  at  various 
points  along  the  a-N  curve.  Values  of  a  at  these  points  were 
used  to  calculate  AK,  the  stress  intensity  factor  range  by  using 
methods  described  earlier  in  pages  6-8.  Values  of  are  plotted  against 

AK  in  figures  43,44  &  45  for  initial  disbond  sizes  o£  19.05  mm 
( 0 . 7 5 K ) ,  25.4  mm  (1 " )  and  31.75  mm  (1.25  in)  respectively.  These 
dimensions  are  those  of  the  diagonals  of  the  diamonu  shaped  dis¬ 
bonds  i.e.  approximate  measure  of  disbonds  at  locations  3  and  6. 

A  least  square  fit  to  the  data  points  in  figures  43,  44  and  45 
yield  the  following  values  of  the  parameters  C  and  n  for  the 

Cj.  y. 

power  law  relationship  ^  *  C(AK) 

C  -  2.044  x  10-23  in  SI  units  (N,m) 

»  2.222  x  10~13  in  lb.  and  inch  units 
n  «  2.776 

It  is  noted  that  the  disbonds  did  not  propagate  in  a  self 
similar  fashion  and  pointwise  measurements  were  used  to  obtain 
the  empirical  crack  propagation  law.  However,  the  relationship 
is  obtained  from  measurements  made  near  the  free  edges,  which 
are  purely  in  a  state  of  plane  stress  and  therefore,  it  is  ex¬ 
pected  that  this  relationship  will  be  useful  for  predicting 
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growth  of  flaws  under  such  a  state  of  stress.  Care  must  be 
taken  in  using  the  data  presented  here  for  studying  propagation 
of  disbonds  under  states  of  plane  strain  or  other  general  con¬ 
ditions  of  loading  in  the  case  of  disbonds  having  two  dimensional 
planar  form. 

BUCKLING  OF  DISBONDED  COMPRESSION  SKIN 

Properties  of  beam  elements  ard  the  honeycomb  core  used  for 
correlations  with  experimental  results  are  shown  in  Table  10. 

Results  of  failure  stress  versus  disbond  length  are  shown  in 
figure  46.  The  analytical  prediction  is  similar  to  the  Euler 
stress  versus  column  slenderness  ratio  curves.  The  cut-off  is 
set  at  the  test  data  for  the  compressive  strength  without  de¬ 
fects,  which  is  equal  to  689.5  MPa.  It  is  clear  that  the  analyti¬ 
cal  predictions  correlate  extremely  well  with  experimental  data 
when  the  disbond  length  is  19.05  mm  or  more. 

Experimental  results  for  12.7  mm  long  disbonds  (about  414  MPa) , 
however,  fall  well  below  the  compressive  strength  of  689.5  MPa 
and  the  buckling  stress  of  793  MPa.  This  discrepancy  appears  to 
be  due  to  inelastic  behavior.  Ttst  data  for  0°  layers  show  that 
the  axial  stress-strain  response  is  linear  up  to  the  strain  level 
under  consideration.  +45°  laminates  show  a  slight  deviation  from 
linearity,  but  it  is  likely  to  have  no  influence  on  the  axial 
and  bending  stiffnesses  of  the  beam  elements,  which  are  dominated 
by  the  0®  layers.  It  is  known,  however,  that  shear  stress-strain 
responses  of  fiber  reinforced  composites  are  highly  nonlinear 
and  therefore,  in  the  presence  of  axial  compressive  stress,  the 
value  of  tangent  modulus  in  shear  (ratio  of  shear  stress  incre¬ 
ment  or  strain  increment)  could  be  significantly  less  them  the 
initial  modulus  (especially  in  +45°  layers).  Estimation  of  the 
correct  value  of  tangent  modulus  will  involve  the  formulation 
of  an  appropriate  yield  criterion  and  a  trial  and  error  procedure 
for  determining  the  value  of  critical  stress.  Such  estimation 
is  beyond  the  scope  of  this  work.  To  study  the  effect  of  reduced 
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shear  modulus,  calculations  were  repeated  (for  12.7  mm  disbond 
length)  with  a  shear  modulus  1/4 th  of  the  initial  modulus.  The 
results  show  that  the  buckling  stress  is  reduced  to  541  MPa  from 
793  MPa.  Lower  values  of  tangent  modulus  will  yield  critical 
stress  comparable  to  414  MPa  observed  experimentally.  It  appears, 
therefore,  that  the  analytical  flaw  criticality  criterion  is  in 
good  agreement  with  experimental  results  for  elastic  material 
behavior. 
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DISCUSSIONS  AND  CONCLUSIONS 

In  the  analytics 1/experimental  data  correlation  studies 
presented  herein,  usefulness  of  the  proposed  NDE  methodology 
has  been  demonstrated.  The  following  important  conclusions 
and  comments  can  be  made  based  on  the  results  of  the  present 
study. 

1.  Use  of  ultrasonic  "C"  scans  is  an  effective  technique 
for  detection  and  growth  of  the  types  of  interlaminar  disbond3 
considered  in  this  study. 

2.  Linear  elastic  fracture  mechanics  approach  in  con¬ 
junction  with  the  methods  of  stress  analysis  presented  here  is 
useful  for  assessing  criticality  of  flaws  in  presence  disbonds 
which  are  in  a  state  of  shear.  The  analysis  method  can  be  ex¬ 
tended  to  the  case  of  multiple  number  of  flaws  at  different 
locations  as  in  the  case  of  impact  damage. 

3.  A  semi-empirical  crack  growth  law  can  be  used  for  pre¬ 
dicting  growth  of  flaws  due  to  cyclic  shear.  It  appears  that  by 

the  use  of  this  law  and  the  failure  models  proposed  herein  an  assess¬ 
ment  of  residual  strength  and/or  residual  lifetime  can  be  made. 
Further  study  in  this  direction  should  help  in  checking  the  ac¬ 
curacy  of  such  assessments. 

4.  Presence  of  disbonds  can  significantly  lower  the  load 
carrying  capacity  of  a  laminated  compression  member  and  reduced 
(or  residual)  strength  can  be  estimated  using  the  model  used 
here. 

5.  It  is  noted  that  all  methods,  results  and  criteria  de¬ 
veloped  herein  are  of  limited  use,  i.e.  for  the  types  of  disbonds 
and  geometry  of  the  beam  type  structure  considered.  Use  of 
these  results  for  other  types  of  disbonds  (2D  planar  form  is  an 
example)  and  different  types  of  structure  should  be  made  with 
due  consideration  to  the  difference  in  such  geometric  config¬ 
uration  and  further  studies  are  needed  before  standard  tools 
for  practical  use  are  developed. 

6.  Moisture  conditioning  has  a  significant  influence  on 
attenuation  of  ultrasonic  waves  and  damages  induced  by  such  con¬ 
ditioning  may  be  assessed  from  the  results  of  wave  propagation 
studies  in  conjunction  with  the  analytical  relationships  presented. 
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Table  1 

Static  Shear  Test  Results 
for  Precracked  [ (0/90) fi012l g  Laminates 


Sample 

No. 

Defect 

-Type 

Precrack  Length 
mm  (inch) 

Failure  Loan 
N  (lb.) 

1-1 

1"  Square 

70  (2.75) 

1451  (326.4) 

1-3 

1"  Square 

30  (1.13) 

427?.  (961.2) 

1-6 

1"  Square 

26  (1.02) 

4233  (952.4) 

D-2 

Diamond 

39  (1.54) 

2940  (661.4) 

D-5 

Diamond 

46  (1.81) 

2332  (524.8) 

D-3 

Diamond 

40  (1.57) 

2204  (496.0) 
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Table  2 

Static  Shear  Test  Results  for  [  (0j|/±452/+452/04)g]  g 
Laminates  with  Diamond  Shaped  Defects 


Specimen 

Debond  Length 
mm  (in.) 

Failure  Load 
nITETI 

Average 

"nTIETT 

0.75-1 

19.05(0.75) 

9021(2028) 

0.75-2 

19.05(0.75) 

9314 (2094) 

9167(2061) 

1.00-1 

25 . 4 (1.00) 

8433(1896) 

1.00-2 

25.4(1.00) 

8611(1936) 

1.00-7 

25.4(1.00) 

8184(1840) 

1.00-8 

25.4(1.00) 

8331(1873) 

8389  (1«*86) 

1.25-1 

31.75(1.25) 

7255(1631) 

1.25-2 

31.75(1.25 

5685(1278) 

1.25-5 

31.75(1.25) 

7842(1763) 

1.25-6 

31.75(1.25) 

7450(1675) 

7059(1587) 
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Table  3 

Fatigue  Test  Results  for  t04/+452/+452/04) g] g 
Laminates  with  Diamond  Shaped  Defects 


Laminate : 


Specimen 

0.75-7 

Debond  Length 

S  Level 

Cycles  to  Failure 

mm  ( in . ) 

19.05(0.75) 

0.6 

35,045 

0.75-8 

19.05(0.75) 

0.6 

12,035 

0.75-3 

19.05(0.75) 

0.5 

68,121 

0.75-4 

19.05(0.75) 

0.5 

64,433 

0.75-5 

19.05 (0.75) 

0.4 

0.75-6 

19.05(0.75) 

0.4 

725,293 

1.00-9 

25.4(1.00) 

0.6 

9,627 

1.00-10 

25.4(1.00) 

0.6 

11,199 

1.00-4 

25.4 (1.00) 

0.5 

46,070 

1.00-6 

25.4(1.00) 

0.5 

79,991 

1.00-11 

25.4(1.00) 

0.5 

48,751 

1.00-12 

25.4(1.00) 

0.5 

48,835 

1.00-13 

25.4(1.00) 

0.4 

150,623 

1.00-14 

25.4(1.00) 

0.4 

194,019 

1.25-3 

31.75(1.25) 

0.75 

5,034 

1.25-4 

31.75(1.25) 

0.75 

5,000 

1.25-7 

31.75(1.25) 

0.5 

27,572 

1.25-8 

31.75(1.25) 

0.5 

29,000 

1.25-9 

31.75(1.25) 

0.5 

28,468 

1.25-10 

31.75(1.25) 

0.5 

16,095 

1.25-11 

31.75(1.25) 

0.4 

47,477 

1.25-12 

31.75(1.25) 

0.4 

57,633 

1.25-13 

31.75(1.25) 

0.3 

>2,250,000 

1.25-14 

31.75(1.25) 

0.3 

412,957 
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Table  4 


'Test  Results  for  Buckling 


Sample 

Defect 

Location 

Number 

Size (in.)  mm 

Location 

18 

None 

19 

None 

— 

20 

None 

— 

21 

None 

— 

IT 

(1/2)  12.7 

Center 

10B 

(1/2)  12.7 

Center 

1ST 

(3/4)  19.05 

Center 

16T 

(3/4)  19.05 

Center 

16B 

(2/4)  19.05 

Center 

17T 

(3/4)  19.05 

Center 

17B 

(3/4)  19.05 

Center 

3T 

(  1  )  25.4 

Center 

7B 

(  1  )  25.4 

Center 

12T 

(  1  )  25.4 

Center 

8B 

(  1  )  25.4 

Center 

6T 

(1.5)  38.1 

Center 

5T 

(1.5)  38.1 

Center 

11T 

(1.5)  38.1 

Center 

9B 

(1.5)  38.1 

Center 

3B 

(  1  )  25.4 

Near  Surface 

12B 

(  1  )  25.4 

Near  Surface 

7T 

(  1  )  25.4 

Near  Surface 

8T 

(  1  )  25.4 

Near  Surface 

5B 

(1.5)  38.1 

Near  Surface 

6B 

(1.5)  38.1 

Near  Surface 

11B 

(1.5)  38.1 

Near  Surface 

9T 

(1.5)  38.1 

Near  Surface 

i 
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0 

3 

1 


Load  at  Stress  at 


Failure,  N 

Failure (MPa) 

9,709 

642 

10,101 

668 

10,983 

726 

11,670 

771 

Avg. 

702 

5,982 

395  j 

6,571 

415 

Avg. 

405 

6,963 

460  i 

4,560 

301 

5,492 

364 

) 

4,707 

311 

7,183 

480 

Avg. 

383 

1 

2,324 

154 

i 

2,942 

194 

2,599 

172 

i 

2,942 

194  j 

Avg. 

179  j  ; 

1,314 

86.9  j  j 

1,373 

90.8  1 

1,304 

86.2 

i 

1,618 

107 

1 

Avg. 

92.7 

i 

98 

6.48 

i 

98 

6.48 

i 

218 

14.3 

i 

i 

128 

8.43  *  j  : 

— 

Avg. 

8.92  j  | 

— 

i 

i 
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Table  5 


Moisture  Conditioning  Effects  on  [±.45]  2g  Laminates 


Specimen 

Conditions 

Shear 

GPa 

Modulus 

(msi) 

Shear 

MPa 

Strength 

(ksi) 

422 

Pristine 

5.6 

(.81) 

82.7 

(12.0) 

421 

Pristine 

5.4 

(.78) 

86.2 

(12.5) 

423 

Pristine 

5.5 

(,l!> 

81.4 

(11.8) 

Avg .  5.5 

(.8) 

Avg.  83.4 

(12.1) 

424 

120  days 

8 

100% 

RH 

4.9 

(.71) 

71.7 

(10.4) 

411 

120  days 

8 

100% 

RH 

4.6 

(.66) 

67.6 

(9.8) 

414 

120  days 

8 

100% 

RH 

4.6 

(.67) 

71.0 

(10.3) 

Avg. 

4.7 

(.68) 

Avg.  70.1 

(10.2) 
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Table  6 


Pulse-Echo  Attenuation  Coefficient  a 


Sample  ID 

Thickness  ram 

a  ram”1 

Dry  Resin  No.  1 

5.08 

.106 

Dry  Resin  No.  2 

6.84 

.110 

Dry  Resin  No.  3 

5.08 

.124 

Dry  Graphite 
±45  Graphite  No. ID 

1.15 

.0190 

+45  Graphite  No. 2D 

1.15 

.00097 

+45  Graphite  No. 3D 

1.15 

.0072 

Saturated  Graphite 

+45  Graphite  No.lW 

1.15 

.058 

+45  Graphite  No.2W 

1.15 

.079 

+45  Graphite  No.3W 

1.15 

.026 

Thick  Beer,  No.  1 

8.86 

.0136 

Thick  No. 2 

8.86 

.0101 

38- 
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Table  7 

Stress  Intensity  Factors  and  Strain  Energy  Release  Rates- 
n04/+452/+452/04)  g]  s  Laminates 


Elasticity 

a/A 

KIc2H 

KJd2H 

ft 

M 

M 

K 

G  EqA 

$ 

4 

Theory 

P  /a 

P  7a 

P  /a 

P* 

$ 

•4 

at  c  &  d 

at  c  &  d 

j 

.004 

-.0005 

.002 

-1.06 

.206 

.025 

-.005 

.006  * 

-1.09 

1.38 

i 

Averaged 

.0417 

-.004 

.005 

-1.15 

2.55 

i 

Properties 

.0833 

.003 

-.004 

-1.30 

6.53 

! 

\ 

i 

.1667 

.0004 

.0001 

-1.60 

19.8 

i 

1 

i 

• 

.26 

.0001 

.004 

-1.89 

42.9 

j 

j 

.004 

-.0004 

.002 

-1.02 

-139 

j 

i ' 

.025 

-.005 

.006 

-1.08 

.976 

; 

f 

Considering 

.0417 

'  -.004 

.005 

-1.17 

1.89 

Stacking 

.0833 

.003 

-.003 

-1.37 

5.20  j 

Sequence 

.1667 

.0004 

.0003 

-1.72 

j 

16*5  | 

i 

j 

.26 

.0002 

.004 

-2.05 

36.4  | 

\ 

! 
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Table  8 

Critical  Loads  and  Strain  Energy  Release  Rates  from 
Various  Theories-I (04/+452/+452/04)sJ  gSamples 


a/e 

hMKESSSIH 

an:TT'~n 

KM 

in 

153 

217 

n 

m 

HKHi 

■uu 

Specimen 

A/l 

Er/2,N 

(Sxptl . ) 

_ Ge»  N/m _ 1 

Elasticity 

(average) 

Elasticity 

(exact) 

Beam 
(Ref .1] 

1-1 

3187 

1185 

986 

864 

1-2 

1  /& 

3383 

1335 

1111 

974 

1-3 

A/  O 

2304 

619 

515 

452 

1-4 

25*8 

787 

655 

574 

Average 

9B2 

TTT 

mm 

i  mi— 

wmiMmm 

WLi+m 

1/2-2 

1/12 

k'H 

785 

624 

432 

1/2-3 

882 

702 

486 

Average 

mam 

TO 

664 

460 

I 


I 
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Table  9 

Critical  Loads  and  Strain  Energy  Release  Rates 
from  Precracked  [ (0/90) fi012] s  Laminates 


Sample 

a Ik 

1-1 

.46 

1-3 

.197 

1-6 

.171 

D-2 

.256 

D-5 

.302 

D-3 

.262 

\ 


Par/2, N  <JC,  N/m 

Elasticity  (Exact) 


726 

765 

2,138 

1,144 

2,118 

879 

1,471 

916 

1,167 

725 

1,103 

513 

Avg.  824 


| 

j 


► 

t 


f 
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Table  10 

Properties  of  Beam  Elements  and  Core 
of  Sandwich  Beams 


Beam  Stiffness 

Beams  1  and  2 

(0/+45/+45/0) 

Beams  3  and  4 
(0/+45/+45/0)# 

Axial  a 

4.22  x  107  N/m 

8.44  x  107  N/m 

Bending  D 

3.444  N.m 

19.14  N.m 

Coupling  3 

0 

0 

* 

Shaar  k55C55 

3.20  x  106  H/m 

6.40  x  106  N/m 

thickness 

.7601  mm 

1.5202  mm 

Core  (Honeycomb  ALC- 1/8-5052-. 006-22.1  ~  Aluminum  Corrugated) 
Properties 
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figure  1.  Laminated  Beam  Containing  Disbonds 
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Figure  5.  Disbond  Fabrication  Method 
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.x-Z'. 
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Figure  18.  Sandwich  Beam  Specimen  Geometry 


*  6 , 963N 
»  460MPa 


1 

4 


6 


e 


c' 


Figure  21.  ^  vs, 
19. 05mn. 


*c  for  Specimen  15T 
Center  Defect 


NADC 


ec,  10'3 


Figure  22.  vs.  eQ  for  Specimen  16T, 
19.05mm.  Center  Defect 


■  2,942N 


1 


e 


c' 


Figure  23.  eT  vs.  ec  for  Specimen  7B 
25.4mm.  Center  Defect 


Figure  24.  eT  vs.  eQ  for  Specimen  5T 
38.1mm.  Center  Defect 


218N 


Figure  26.  eT  vs.  eQ  for  Specimen  9T 


38.1mm.  Near  Surface  Defect 
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Figure  37.  Critical  Load  for  Varying  Disbond  Length 


10  20  30  40  50  60 


N,  1000  Cycles 


Figure  39.  a  vs.  N  for  Specimens  0.75-3  and  0.75-4,  S»0.5 


12  (1) 
11  (5) 
11  (2) 


Curve  Fit 


□  1.00-10 
O  1.00-9 
&  1.00-14 
+  1.00-13 

X  1.00-12 
0  1.00-11 
^  1.00-6 
O  1.00-4 


Figure  46.  Data  Correlation  for  Buckling  Tests 
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APPENDIX  A-l.  STRESS  ANALYSIS  OF  DELAMINATED  BEAM  IN  SHEAR 
SERIES  SOLUTION 


Consider  a  laminated  beam  of  span  i  as  shown  in  Figure  1 
under  a  state  of  generalized  plane  stress  or  plane  strain  such 
that  the  displacements  in  all  layers  can  be  described  by  two 
components  U(«u4)  and  W(-w£),  The  differential  equation  in  terms 
of  the  nondimensionalized  displacements  and  w^  of  layer  i 


are  as  follows: 

.2 


(C 


13 


+  C55> 


3x 

a\ 

3x3z 


ui  + 


(C 


13 


+  c^5) 


+  C 


55 


*2wi 
Sx3  7 


— V"! 

3x* 


(A-l.l) 


,1 
'jk 

plane  stress. 


where  Ct,.  are  the  effective  moduli  of  layer  i  in  the  case  of 


Assume  solutions  of  (A-l.l)  in  the  form: 


09 

u.  *  I  u.  (z)  cos  mrrx 
1  m*l  im 


430 

w.  *  E  w.  (z)  sin  imrx 

i  im 

m=  J. 


(A-l. 2) 


Substituting  (A-l. 2)  in  (A-l.l)  and  solving  for  u.  (z)  and  w.  (z) ,  one 

im  im  .  3.m 

can  express  the  stresses  c  and  t  ,  m  layer  i  corresponding 

7,  Z  ZX 

to  the  roth  harmonic,  i.e.. 


o  *  E 

cr* 

sin 

roirx 

zz 

ZZ 

m 

(A-l. 3) 

T  *  I 

T 

cos 

mrrx 

ZX 

ZX 

in  terms  the  displacement  components  at  interfaces  i  and  i+1. 
In  particular  the  boundary  tractions  at  z«±h/2  %  can  be 
written  as: 


-8*- 
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where  +  and  -  signs  indicate  values  of  o™z  and  at  z=4h/2 
and  -h/2, respectively.  A™fl (a,P-l,2)  are  2x2  stiffness  matri¬ 
ces,  each  element  of  which  is  a  complicated  function  of  the 
layer  stiffnesses  and  its  thickness.  Detailed  expressions 
for  anisotropic  layers  are  given  in  table  A-l-1.  Writing 
the  displacement  vector  twim>uimlT  as  an<*  traction  vector 

— [om,Tm]T  as  Sm,  one  obtains  the  following  equations  from 
IMT 

boundary  and  interface  conditions  for  the  layered  system  in 

terms  of  the  prescribed  tractions  and  S™  on  top  and  bottom 

+  n  — *m 

surfaces  of  the  beam  and  the  displacement  discontinuity  u 
at  the  interface  (i+1)  containing  the  disbonds  [see  equation 
(A-l. 13)3 .  Note  that  superscripts  (+)  and  (-)  indicate  quanti¬ 
ties  on  the  top  and  bottom  of  interfaces  (i+1)  respectively. 


Am(1,u 

An  ux 


m(l)-m 

a,  2  «2 


(A-1.5) 


NADC-77278-30 


NADC-77278--30 


(A-1.17) 


(A-1.18) 


(A-1.19) 


vrhich  satisfy  (7  -1.13)  and  reduce  the  system  (A-l.5-7),  (A-1.15) 
and  (A-l.  10-12)  ^.o: 


91- 
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where  [K  ]  is  the  global  stiffness  matrix  for  the  laminate 
m 

and: 


5‘n(k)  -  Am(k>  -  A(k)*;  a=l,2  and  k=i,  i+1 
oo  oa  oa 

One  should  note  that  [A™^k*  ]  as  well  as  [Amjk^  ]  (a?<$) 

oa  op 

approach  [0]  as  m  is  increased. 


(A-1.21) 
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REDUCTION  TO  A  'SET  OF  INTEGRAL  EQUATIONS 

Now  we  introduce  two  unknown  functions  such  that  the  dis¬ 
placement  discontinuities  at  interface  (i+1)  are  given  by: 


u(+)  -  u(‘-) 

Ui+1<*>  i+l<*> 


;  0  <  x  <  c,  d<x<l/2 


-u 


@2  (x)  dx;  c  _<  x  <_  d 


(+1  <-l  n 

w  (x)  “  w  (x)  “  0 
i+lw  i+1  'x' 


•  ■  */: 


;  0  <  x  <  c,  d<x<l/2 


0^(x)dx;  c  <_  x  <,  d  (A-l 


such  that: 

d 
’  c 


(x)dx  =  0;  i=l,2 


(A-l 


and 


u(l-x)  =  -  u(x)  ;  0  <_  x  <_  1/2 
and  w(l-x)  =  w(x) 


(A-l 


Equation  (A-l. 23}  applies  to  a  structure  loaded  symmetrically 
about  the  center  line  x«l/2.  Using  (A-l. 23)  and  (A-l. 22)  u  in 
equation  (A-l. 13)  can  be  expressed  as: 


r  ^ 

’  w*m 

1 

SB  - 

/  9-.(y)  cos  miry  dy 

-'c  A 

1 

e 

* 

a 

_ 1 

mn 

f 

/  6-<y)  sin  miry  dy 

(A-l 


.22a) 


.22b) 


.23) 


.24) 
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Solving  for  displacement  components  (j*l, n)  ,  substituting 
their  values  as  well  as  (A-1.19)  and  (A-1.24)  in  (A-1.8)  and  sum¬ 
ming  the  resulting  expressions  for  from  m-1,3.. »,  one  can 

express  the  stress  free  conditions  on  top  surfaces  of  the  disbond 
in  terms  of  a  set  of  coupled  pair  of  integral  equations. 


i(d  v 
1  y‘ 


y)  H 


12 


x  ay  -  Sr  Vx)  V 

Hu 


(x»y)6j (y)dy  *  PL(x) 


(A-l. 25a) 


1  /  e2(y>  H12  2  fA  • 

w  -rnr  d»  +  #  Vx>  +2I 

yc  22  3  ±JC 

00 

•  11  (  -s/2 

*  -  —  [^—  +  2J  e  '  sech  s/2  cosh  sy  sinh  sx  ds] 


~  Z  sin  mirx  cos  miry 


H 


11 


*  4  in 

K. -  -  — — -  Z  H.  _  sm  mirx  sxn  miry 
12  „o  12 

H11 


*  4  in 

K-,  *  — - —  Z  H_.  cos  mirx  cos  miry 
21  o  21  .  J 

H22 


"22 


i  (4-+  2  f 

it  1  x+y  J 


-3/2 


e  sech  s/2  sinh  sy  cosh  sx  ds] 
o  (A-*l.  25b) 

+  Z  11^2  cos  mirx  sin  miry 


H 


22 


pi  *  "Tr  2  Plm  sin  mirx 
x  „o  im 

Hn 


P  4  *.  avin 

2  2  P22  008  m1TX 
H22 
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The  constants  H°x,  H°2,  H°2#  H™g(a,8-1,2)  and  P  and  ?2m 
are  known  mid  the  summation  sign  indicates  sum  over  m»l,3,...«. 
One  should  note  that  use  has  been  made  of  some  well  known  identi¬ 
ties  (ref.  18)  to  obtain  Equations  (A-1.25a,b).  In  (A-1.25) 

r 

Kij(x,y)  are  the  regular  kernels  and  Cauchy  kernels  appear  in 
the  first  terms.  These  types  of  equations  have  been  considered 
in  reference  19.  By  the  substitutions  given  below: 

x  “  a1t  +  a2 
y  “  a^t  +  »2 
ax  «  (d~c)/2 
a2  *  (d+c)/2 

0i(T)  =  A  (fc)  *  P^x) 


Equations  (A-1.25)  can  be  reduced  to: 


i/1  -■ 

*  -i 


6^(1)  ^12  ^ 

dx - ^  6„(t)  +  z 


lll 


j’ 


r1 

1  J-i 


1  I1  <Mt)  K12 

r-  /  -2—-~  dx - H  6.(t)  + 

IT  A  t-t  „0  1 


-1 


H 


22 


2  rl 

1  J 

3-1 


KXj  (t,x)6..  (x)dc 


K2j (t,T) 0^ (T)dT 


px(t) 

(A-1.26) 

P2(t) 


Following  methods  given  in  reference  20,  one  can  study  the 
singular  nature  of  functions  0^  and  02  as  t  +  ±1.  When  the  pro¬ 
perties  of  layers  i  and  (i+1)  are  identical  H°2-0  and  consid¬ 
eration  of  Cuchy  kernels  yields  the  well  known  square  root  type 
of  singularity.  When  H°2  ^  0  the  index  of  the  singularity  is 
of  the  type  -  l/2iit ,  which  is  the  characteristic  stress  singu¬ 
larity  at  the  tip  of  a  cracx  between  two  dissimilar  materials. 

The  system  of  Equations  (A-1.26)  and  (A-1.22b)  must  be  solved 
numerically  and  the  suitable  methods  of  solution  are  different 
for  the  cases  (i)  HX2»0  and  (ii)  In  this  study  we  will 
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consider  H12-0  and  employ  the  colocation  method  outlined  in  ref¬ 
erence  19  to  reduce  the  equations  to  a  system  of  algebraic  equa¬ 
tions  for  determining  the  functions  ^(t),  1*1,2  defined  by: 

Qi(x)  -  (1-tV1/?  ♦i(T).  1*1,2  (A-l.,27) 


By  solving  the  set  of  equations  for  $.(a)  at  points 

7T  f  2k“ }  )  ^ 

i^cos  — —  ;  (k*l,2...n)  can  be  evaluated.  The  number  n 
denotes  the  number  of  collocation  points  used  for  discretizing 
the  system.  The  stress  intensity  factors  at  x*c  and  d  corres¬ 
ponding  mode  I  and  II  are  given  by: 


KId  "■ 

•H° 

H11 

Vd^ 

h 

(l)/4  sj~2 

K!c  " 

H?i 

v'd^c 

♦l 

(-l)/4  s [2 

KUd  - 

H 22 

J  d-c 

^2 

<l)/4  sJT 

KIIc 

-H22 

Jd-c 

^2 

(-D/4  VT 

(A-1.28) 


and  the  strain  energy  release  rate  can  be  shown  to  be: 


G 


KII  1 

„o  J 


(A-1.29) 
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TABLE  A- 1-1 

STIFFNESS  MATRIX  FOR  A  LAYER 

Let  A  -  C55  C33 

B  "  (C13  +  C55)2  “  C11  C33  "  C55 

C  "  C11  C55 
D  -  B2  -  4AC 

For  the  type  of  anisotropic  layers  considered  B<0  and  D>0. 
The  stiffness  matrix  for  harmonic  m, 

’An  A12- 
tKl  “  A  A 

La21  a22J 

is  given  1  7: 

-  2 (b, -b_) (a,  coth  p  ct,h  -  a,  coth  p  ci„h)/A* 

11  121  ml  2  m2m 

k12  “  t2(dIal+d2a2)‘2(d3.a2+d2al,^m3/Am 

k.,  =  2(b.-b_) (-a,  cosech  p  a,h+a„  cosech  p  a_h) / A* 

13  12  l  m  1  2  m2  m 

k14  *  2  (d2a1-d1a2)  j-(tanh  p^a^h/j  coth  pma2h/2~  cbth 
pmalh/2  tanh  pma2h/2^ 

k22  -  2(d2a1-d1a2)(ax  coth  p^h-aj  coth  P^W/A* 

k24  -  2*d2al“dla2*  ^-al  cosech  pma2h+a2  ccsech  Praalh3//A 
k33  *  kll 
k2l  “  k12 
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TABLE  A-l-1  (Continued) 

k34  "  k43  -  "k12 

k41  "  k14 

k23  "  k32  "  "k14 

k31  "  k13 

k42  "  k24 

k44  “  k22 

Where  p  ■  mir 
in 

rm  *  J  [tanhlp^jh/j)  ooth  IPm»2fl/2l+  Coth(pma1h/21  tanh 

<pm“2h/2>] 


4S  -  2<»12  +  *2  - 

ax  =  { (-B  +  JD)/2A}1/2 

<*2  -  (  (~B  -  s'd)/2a}1/2 

»i  *  '  +  CU>'<C13  +  C55>“i 

bi  ’  C33*i“i  -  C13 

di  '  C55<ai  +  V 

ei  '  C13Vl  -  CU 


If  D<_0  different  expressions  for  the  elements  k 


ij 


hold. 


I 
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APPENDIX  A- 2.  FORMULATION  OF  THE  BUCKLING  PROBLEM 

For  the  Beams  1,  2  and  3  (i*l,2,3)  define  the  following 
nondimens ional  quantities. 


Coordinate 
Thickness  H 


(i) 

(i) 


X(i)/i 

h(i)/2! 


Midplane  displacement  in  x-direction  U 


(i) 


Beam  deflection  in  z-direction 

Bending  rotation 

Bending  stiffness 

Axial  extensional  stiffness 

Coupling  stiffness 
( ex tens ion-bending ) 

Shear  stiffness  (effective) 
Foundation  modulus  in  shear 
Foundation  modulus  in  extension 
Axial  Force 


W 


(i) 


w(i)A 


*(i) 

d^> 

a(i) 

b(i) 

k(i) 

s(i> 

tU) 


»  D 
*  A 
=  B 

■  k 

=  s 


,<i> 

11 

/0<i> 

(i) 

11 

(i) 

11 

55C 

ggi) i2/D 

(i) 

•H’ 

,(i) 

(A-2.1) 


11 


8(i)n^S(i)<(1)i2/D^ 


*  (i) 

Nx  is  the  axial  force  in  beam  1.  For  beam  1  foundation  moduli 
are  zero.  Governing  differential  equations  in  nondimensional  form 
for  the  laminated  beams  are  obtained  from  laminated  plate  theory 
including  effects  of  shear  deformation,  i.e.: 


a(i)U(ti  +  b(l)il>(**  -  sU)  (i)H(l)  )  =  0 


,xx 


.XX 


b^U*i)  +  d^ij;^  +  W^)  +  s (U  =  0 


,xx 


,xx 


k(i)  -  t(i)W(i)  -  3(i)n*W  -  0 

•  *  »xx  ,xx 


(A-2.2) 
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Stiffness  Matrix  for  Beam  1 

For  beam  1,  *  0  and  the  nonzero  characteristic 

roots  of  (A-2.2)  are  given  by: 


A  -  ±iA0 


where  AQ 


ak8n 

T, 


(ad-b  ) (k-0n  ) 


(A-2.3) 


Superscript  (1)  has  been  omitted  in  equation  (A-2.3)  and  such  super 
scripts  will  be  omitted  except  where  necessary.  Considering  the 
deformation  patterns  which  are  symmetric  about  ■  0, 

the  nondimensionalized  stress  resultants  at  x^^  »  1/2  are 
related  to  the  displacements  at  that  point  in  terms  of  the  stiff¬ 
ness  matrix  i.e. 


'nx 

u 

mx 

-  IK(1>] 

qx 

x=l/2 

w 

x-1/2 


(A-2.4) 


where  n  ,  m  are  the  axial  force  and  bending  moment,  respectively, 
/  *  *  *  [1  ) 
and  qx  ^  qx  ”  3x'  qx  bein9  the  shear  force.  K  7  is  given  by: 


[K(1)] 


2a 

2b 

0 


2b2/a  +  (d-b2/a) Aq  cot  AQ/2  0 


(A-2.5) 


-100- 


NADC-7 72 78-30 


Characteristic  Equations  for  Beams  2  and  3 

Characteristic  roots  X  for  equation  (A-2.2)  are  the  solutions 

of: 

X6[(k-8n*) (ad-b2)] 

+  X4  [ak2-t  (ad-b2)  -  (k-8n*)  {a(k-*H2)  +  sd  4  2beH}]  (A-2.6) 

+  X2 [t{a (k+sH2)  +  sd  +  2bsH}  -  8n*sk]  -  skt  *  0 


Once  again  the  superscript  i  is  omitted  in  equation  (A' 2. 6)  for 

2 

beams  with  geometry  and  stiffnesses  under  consideration  X  has 
three  solutions,  i.e.  : 

2 

(i)  one  real  positive  root,  X^ 

(ii)  two  complex  roots,  one  conjugate  of  the  other, 

X2  3  «  (o2-a2)  ±2i  (A-2.7) 


Stiffness  Matrix  for  Beam  2 

For  beam  2  we  consider  deformation  patterns  symmetric  about 
0,  and  introducing  three  integration  constants,  E^E^jE-j  (*E) , 
the  displacements  [U]  *  [U  ,ip  ,w]  x(2)  _  and  stress  resultants 
[N]  -  x<2)  m  2_/2  are  expressed  as  : 


tU]  -  [F ]  [E] 
[N]  =  [G]  [E] 


(A-2.8) 


Matrices  F  and  G  are  given  by  : 
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nQ  cosh  X^/2  ^nlrl“n2r2^  *n2rl+nlr2^ 

G  ■  mQ  cosh  X^/2  *mlrl“K2r2*  *in2rl",’mlr2* 

qQ  sinh  X^/2  (qlpl"q2?2)  (q2pl+qlp2) 

(A-2.10) 

2 

where  X1  is  the  positive  root  of  X^  and  and  a2  are  positive. 

Also,  =  - (bX?+sH) / (aX2-s) 

U  i  X 

f 0  =  -k/{ (k-8n*) X2  -  t} 
nQ  -  (aeQ+b)X1 

m0  -  (b«Vd)Xl 

qQ  =  tf 

Gl  *  -[{b(o2-a2)  +  sH}U<a2-a2)  -  s)  +  4ab  ajoi2]/^ 


e?  =  20^28(5+201) /A  ^ 

=  [a(a2-a2)-  si  2  +  4a2  ct2a2 
fx  =  -k [ (k-Bn  ) (a2+a2)  -  +  ]/A2 
f2  *  k[(k-Sn*) (a2+a2)  +  t]/A2 
A2  *  [{ <k-0n*) <a2-a2)  -  t}2  +  4(k-0n*)2  a2a2] 
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P1  -  sinh  a^/2  cos  a2/2 
P2  =  cosh  a±/2  sin  a2/2 

=  cosh  a^./2  cos  o2/2  ’ 

r2  =  sinh  a^/2  ,in  a2/2 
n^  *  a(e1a1-e2a2)  +  ba^ 

(A-2.11) 

n2  «  a(e2o1+e1o2)  +  bo2 
m^  =  b {e^a^-e2o2)  +  do1 

m2  =  b  (eja^+e^j)  +  da2 

*  tk[(k-0n  )  (a2-a2)  +  t]/A2 
q2  *  2(k-f!n  )  tk/A2 

The  stiffness  matrix  for  beam  2  relating  stress  resultants  and 
displacements  at  x2  =  1/2,  i.e.s 

[N]  «  [K(2)][‘J]  (A-2.12) 

is  therefore  given  by  : 

[K(2)]  =  £GJtF]_1  (A-2.13) 

Stiffness  Matrix  for  Beam  3 

Reversing  the  signs  of  stress  resultants  at  x^  =  0  (to  con¬ 
form  to  the  sign  convention  at  right  hand  end  of  beams  1  and  2) 
and  introducing  three  unknown  constants,  the  displacements  and 
stresses  ?re  given  by  (A-2.8)  where: 
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1  , 


0X1 


fl“l 


f2°2 


(A-2.14) 


(A-2.15) 


2 

Where  is  the  negative  root  of  X^  a2  is  positive  and  ct^  is 
negative.  eQ,  fQ/  nQ/  raQ  and  qQ,  e2»  fj,  n1#  n2,  n^, 
m2  sure  the  same  as  in  (A-2.11)  . 

The  stiffness  matrix  is  given  by: 

[K (3) j  -  [G] (F]"1 *  (A-2.16) 

Using  the  compatibility  conditions  : 

U*1*  -  u(3)  +  i)>*3* 

u(2)  -  u(3)  -  <H(3>-H<2))  *<3> 

(A-2. 17) 

*(1)  »  *<*>  »  *»> 

W(1)  «  w(2)  -  w(3) 
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and  the  equilibrium  equations: 


_  ( i)  .  _  (2)  . 
n  +  n.  + 

„<3> 

X  X 

X 

m(i>  +  m<2)  + 

m(3) 

m 

X  X 

X 

'  (1/  A  '  (2) 

• 

q  +  q 

+  <3 
^x 

*  n'1’ 


(3) 


♦  n<21(H,2)-H(3>) 


(A-2.18) 


one  can  obtain: 

[!C]  [U  (3)]  =  0 


.  _  ,  * 
where  K  is  the  global  stiffness  matrix.  Critical  values  of  n 

(or  are  determined  by  equating  the  global  stiffness-matrix 

to  zero. 


3 

i 

4 

a 

i 

I 

? 

i 
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APPENDIX  A-3.  ELASTIC  WAVE  PROPAGATION  IN  A  LAMINATE 

A  typical  laminate  is  shown  in  fig.  47.  The  laminate  is 
referred  to  a  fixed  coordinate  system  x-,,x.,x.,  while  a  typical 
lamina  with  fiber  reinforcement  angle  0  is  referred  to  x^,x2,x3- 
Consider  a  transverse  wave  in  x3  direction  characterized  by 
the  displacements: 

*  u.(x3,t) 

u2  =  u2(x3,t)  (A-3 . 1) 

u3  *  u3(x3'fc) 

It  follows  that  the  nonvanishing  strains  are: 

2e13  *  ul,3 

2e23  =  U2,3  (A-3. 2) 

e33  *  U3 , 3 

In  order  to  establish  the  equations  of  motion  associated  with 
(A-3.1)  it  is  necessary  first  to  find  the  stresses  associated  with 
(A-3. 2).  For  this  purpose  the  stress-strain  relation  of  a  lamina 

t  I  9 

is  first  written  in  its  material  system  x^,x2#x3>  Thus: 
all  "  Cllll  £11  +  CU22  £22  +  C1122  6  33 
°22  =  C1122  ell  +  C2222  e22  +  C2233  e33 
°33  “  C11222  ell  +  C2233  *22  +  °3333  %3 
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°12  "  2CUl2  C12 


I  I  f 

a23  "  2C2323  e23 


I  I  t 

*C2222"C2233)e23 


i  %  • 

c13  “  2C1212  el3 


{A-3.3) 


The  strains  in  (A-3.2)  and  the  strains  in  (A-3.3) 
are  related  by  tensor  transformation  as  defined  by  the  relation 

I  I 

between  the  x^,x2  and  x^/Xj  axes.  It  follows  that: 


t  I  I 


*11 

“  e22 

-  *«  - 0 

e  33 

"  C33 

G13 

*  *13 

cos0  -  e23 

sin6 

e23 

=  e13 

sin0  +  e23 

COS0 

(A-3.4) 


The  surviving  elastic  moduli  in  (A-3.3)  are  thus  : 

C1122  “  1 
C2233  =  k  “  GT 

C3333  k  +  GT  (A-3.5) 

« 

C2323  *  GT 
C1212  *  GA 
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where  k  is  transverse  bilk  modulus  and  G_  and  G,  are  transverse 

T  A 

and  axial  shear  moduli,  respectively. 

Further  tensor  transformation  on  stress  yield  relations  of 
the  9  rotated  lamina  with  respect  to  the  x^,x2  axes  : 


a13 

2C1313  e13 

+  2C1323  e23 

a23 

*  2C1323  e13 

+  2C2323  G 23 

(A-3.6) 

0  33 

=  C3333  e 33 

where 

2  .  2 
G,  cos  6  +  G_  sin  9 
A  T 

2  2 
G.  sm  9  +  G_  cos  0 
A  T 

(A-3.7) 

(Ga~Gt)  cos 8  sin0 
*  +  gt 


C1313  “ 
C2323  “ 
C1323  " 
C3333  ” 


The  only  surviving  pa  :ts  of  .'.:he  equations  of  motion  are  : 


13,3 


PU! 


°23,/3  ’  pU 


O  «.  v 


1 

2 
pU3 


(A-3.8) 


/ 

where  p  is  the  density  and  a  dot  denotes  time  derivative. 
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Insertion  of  (A -3. 2)  into  (A- 3.8)  and  the  resulting  ex¬ 
pressions  into  (A-3.0)  yields  the  following  wave  equations: 


C1313 

U1 , 33  +  C1323 

U2 , 33 

p 

a. 

n 

(a) 

C1323 

U.l,  33  +  C2323 

U2 , 33 

”  ou2 

(b) 

C3333 

U3, 33 

=  P«3 

(c) 

In  addition  to  these,  displacements  and  tractions  must  be  continuous 
at  lamina  interfaces ►  It  is  seen  that  the  displacements  (A-3.1) 
satisfy  continuity.  The  traction  components  at  an  interface  are 

ai3'  023  and  a33  and  are  9f;*-ven  (A-3. 6). 

In  view  of  (A-3. 7)  it  is  seen  that  o33  associated  with  the 
present  solution  is  always  continuous  while  and  a23  are  not 
continuous  for  adjacent  laminae  with  different  6.  It  is  also 
seen  that  the  discontinuity  is  entirely  due  to  difference  be¬ 
tween  the  values  of  G,.  and  G_.  If  this  difference  is  small  as  is 

A  T 

frequently  found  in  practice,  then  the  discontinuity  may  not  ser¬ 
iously  detract  from  the  present  simple  treatment. 

It  is  also  seen  that  the  discontinuity  only  concerns  the 
shear  wave  equations  (A-3.9;a,b) .  It  may  be  argued  that  if  a 
transverse  pulse  is  transmitted  through  a  small  area  of  a  lamin¬ 
ate  then  u2  and  u^  in  the  laminate  plane  may  be  disregarded, 
since  the  wave  speed  in  transverse  x3  direction  is  so  much 
larger,  thus  the  shear  waves  having  insufficient  time  to  develop 
by  the  time  the  transverse  pulse  has  been  reflected  in  opposite 
direction . 

As  a  consequence  of  all  this,  equation  {  A-3. 9)  is  reliable. 

The  wave  speed  associated  with  it  is  : 


P  c 


vw 


k  +  G„ 


(A-3. 10) 


\ 

•It 

i 


i 

i 

i 


j 
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